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O ■ Introduction 
<N 

We extend the well known Holstein-Primakoff method to the case 
of s/(2,M) algebra. The original method is based on investigation of 
C<") ' quantum systems described by the elements of the enveloping algebra 

of Heisenberg-Weyl algebra satisfying commutation relations of sw(2), 
(30 ! see HIP 10; , and is most effective in the study of shell model of the 

nucleus. Various examples of that can be found in |KM91j . Our gener- 
alization is used to the integration of the quantum-optical models such 
as for example the model of the two-level atom coupled to single mode 
radiation, see e.g. [SB84, Suk89j. 

In this paper we combine the s/(2,M) symmetry and the theory of 
the orthogonal polynomials in order to integrate the wide family of 
multiboson Hamiltonians describing the interaction of bosons with non- 
linear medium as well as the boson-boson interaction. Hence among 
other one can use them for modelling the non-linear phenomena in 
^ ■ quantum optics including Kerr-like effect and parametric generation, 

see (EMEnOToH- 

The paper consists of three sections. In the Section [T] we distinguish 
the two-parameter family of one-mode multiboson Hamiltonians (jl.38|) 
defined as a linear combination of the generators of Lie algebra sl(2, M.) 
obtained by solving the difference equations ()1.4|) - (jl.6J) . We define the 
group of the Bogoliubov-like transformations and using it we distin- 
guish the classes of Hamiltonians related to the appropriate classes of 
orthogonal polynomials, see also Appendix [X] At the end of the sec- 
tion, we construct the system of coherent states which allows us to find 
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new, as we suppose, realization of the discrete series representation of 
the group SL(2,R). 

Section El is devoted to solving the class of two-mode Hamiltonians 
(|2.9|) . which are defined by the combination of Casimirs for sl(2,~fil) 
algebras. Similarly to the one-mode case the integration of (J2.9|) is 
carried out by the orthogonal polynomials, i.e by dual Hahn and con- 
tinuous dual Hahn polynomials. 

In the last section we discuss the possible physical interpretation 
of the quantum models governed by the Hamiltonians investigated in 
previous sections. 

1. One-mode systems 
1.1. Multiboson representations of sl(2, M) in bosonic Fock space. 

Let a, a* be the standard annihilation and creation operators satis- 
fying CCR relations and acting in the Hilbert space TL. Let {|ro)}^l 
be the orthonormal basis in 7i consisting of eigenvectors of number 
occupation operator n := a* a. 

For any fixed I G N, we define the multiboson representation of 
s/(2,M) as the triple of operators 

A :=a (n), A := a_(n) a 1 , A+ := (a*)'a_(n) (1.1) 

defined on dense subset T>q consisting of finite combinations of \n) and 
satisfying sZ(2,R) commutation relations: 

[A_,A+] = A , [A , A±] = ±2A ± (1.2) 

with symmetricity conditions: 

A C A*, A c A% A+ c A* . (1.3) 

We denote by A the operator Lie algebra generated by Ao, A_, A + . 

Relations ()1.2|) and ()1.3|) imply that the functions ao, a_ are real- 
valued and satisfy the following difference equations 

(a (n) - a (n - I) - 2)a_(ra - I) = for n > I, (1.4) 

(n + 1); a_(n) — (n — I + 1)/ a 2 _{n — I) = a (n) for n > I, (1.5) 

(n + l)i a 2 _{n) = a (n) for > n < (21.6) 

where (n)i = n(n + 1) . . . {n + 1 — 1) and one has applied the identities 

(a*)^ = (n - I + 1), a*(a*y = (n + 1),. (1.7) 

Without loss of generality we can assume that a_(n) > for n = 
0, 1, . . . (by multiplying basis vectors by ±1 and by removing common 
kernel of A , A_, A + ). 



The solution to the system of difference equations (jl.4|) . (jl.5j) . (jl.fij) 
is of the form 



a {n) 



+ a (n mod I), 



a_(n) 



n + 1 



+ a (^ mod /) 



+ 1 , 



1.8a) 



;i.8b) 



where [x] is the integer part of x. The values ao(r) for r = 0, — 1 
are arbitrary positive constants corresponding to initial conditions on 
the solution a of the difference equation (jl.4j) . 

In order to express the operators A , A_, A + explicitly in terms 
of the creation and annihilation operators let us define the bounded 
operator 



R: = 



I - 1 



i-i 



T expt-— , v 
2 ' ^ exp(^P) - 1 v ' ' 

m=l I 

for / > 1 and R := for 1 = 1, see |H(X)T 04j. This operator acts on 
elements of the basis by 

R \n) = n mod / \n) (1-10) 
and commutes with operators Aq, A_, A + . Thus one has 



1-Kim 
I 
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y(n-R)|n) 



\n) . 



1.11) 



Finally the multiboson representation of s/(2,IR) is given in terms of a 
and a* by 

A = y(n-R) + « (R), (1.12a) 



J_ Q( n - R) + a (R)) Q(n - R) + l) a'. (1.12b) 

a is these formulae is treated as an arbitrary positive function on spec- 
trum of R, i.e. it is defined by initial conditions a (r), r e spec(R) = 
{0,1,.../- 1}, see ffTSj) . 

The formulae ()1.12j) show that the Hilbert space 7i splits 



i-i 



;i.i3) 



r=0 



onto invariant subspaces 

Tt r := span{ \k) r :=\kl + r) \ k E N U {0}}, (1.14) 
which are eigenspaces of R corresponding to the eigenvalue r. 
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Let us observe that \k) r are eigenvectors of A 

A |A;) r = (2A: + ao(r)) \k) r (1.15) 
and A_, A + act on \k) r as weighted shift operators: 

A \k) r = ^fk{k + a (r) - 1) \k- l) r , (1.16) 

A + |fc) r = v/(A; + ao(r))(A; + l) |* + l) p . (1.17) 

So, ()1.13j) gives the decomposition of the multiboson representation of 
s/(2,M) onto irreducible components. 

Let us now present two examples of the multiboson representation 
of s/(2,R). 

Example 1.1. Let us begin by the simplest case 1 = 1. Then R = 
and cco(R-) = c. Therefore 

A = 2n + c, A_ = Vn + ca. (1.18) 

It is s/(2,R) analogue of Holstein-Primakoff mapping for SU(2), see 
|KM91| . 

□ 

Example 1.2. Let us put / = 2 and the special choice of a^, namely 
cco(R.) = \ + R. Then formulae ()1.12|) simplify to: 

A = n+i K - = \* 2 - ( L19 ) 

This representation is associated with the problem of second-harmonic 
generation, see jWT721 IJD92] . 

□ 

The Casimir operator for A is of the form 

CU := \K ~ (A-A+ + A + A_) (1.20) 

and it can be explicitly expressed by 

C^ = ^a (R)(a (R)-2). (1.21) 

Thus Gj\n r = | a o( r ) ( a o( r ) — 2) and each 7i r is contained in some 
eigenspace of C4. For < a.o{r) < 2 or «o( r ) = 2,3,... repre- 
sentation (|1.15j) - ([1.17)) corresponds respectively to the complementary 
and discrete series of unitary representations of the group <SX(2,R), 
see jVK91j . 
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From (|1.15j) it follows that A is closable with closure defined on 

v i-={^2 v n\n)eH ^2n 2 \v n \ 2 < oo> . (1.22) 



. n=0 



n=0 



Equations (|1.16|) - (jl.l7j) imply that A_ and A + are superpositions of 
diagonal operator and shift operator. Both of these diagonal opera- 
tors are closable with closure defined on T>\. Since shift operator is 
bounded we conclude that the closures of all generators of A are de- 
fined on common domain T>\. From now on Ao, A_ and A + will denote 
these closures. Moreover symmetricity conditions (|1.3|) are replaced by 
stronger conditions 

A = A*, A_ = A* + , A + = A*_, (1.23) 

i.e. operator Ao is self-adjoint and A_, A + are mutually adjoint. 

1.2. Bogoliubov-like transformations. 

Let us consider a linear transformation of A preserving commutation 
relations (jl.2j) and conditions (|1.23|) . We will call such transformation 
a Bogoliubov transformation for the algebra s/(2,R) in the analogy to 
Bogoliubov transformations for Heisenberg-Weyl algebra, see |Bog82 



The group of all transformations of this type is isomorphic to the 
group 03 := R x x Z 2 , where R x := R \ {0}, Z 2 = {-1, 1}, with the 
group operation defined by 

(a, a) ■ (d,t) := (ab r7 ,ar). (1.24) 

One shows that the transformation given by an element (a, a) G 03 acts 
on the generators of A in the following way 

b a , CT (A ) := ii^A + ai^(A_+A + ), (1.25) 
1 a\ (l-a)\ (l + a) 2 A 



ba,.(A_) := — A + a- r— ^-A 



a- 



4a 4a 4a 

KAA + ) = ^ r A + aL_^A + + aL__^A_. 

Let us notice that the domain of operators (|1.25J) is also T>\ and the 
Casimir operator ()1.20j) is invariant with respect to the action (|1.25j) . 

There exists the unitary projective representation (a, a) i — > U ajCT in 
H r of the subgroup R + x Z 2 C 03 which implements the action of b a>(T , 
i.e. 

b a , CT (X) = U a , CT XU: i(T (1.26) 

for all X G A on V x . 



In order to prove this we give the explicit formula for U a>a . First of 
all let us observe that basis vectors \k) r in TC r are eigenvectors of A . 
We will show that b ajCT (A ) has the same spectrum as A and we will 
compute its eigenvectors. To this end we observe that 

b aj(T (A ) \k) r = \k - l) r + a k \k) r + b k \k + l) r , (1.27) 

where 

a k = +a (2k + a (r)) (1.28a) 



-^{k + a Q {r)){k+l). (1.28b) 

If a ^ 1 then the formula (|1.27|) is directly related to three term 
recurrence relation 

xP k (x) = P fc -i(2) + a fe P fc (x) + b k P k+1 (x) (1.29) 

which is valid for any orthonormal polynomials family {-Pn}^L for apro- 
priate choice of coefficients. 

In case when coefficients a k and b k are given by (jl.28j) we obtain 
Meixner polynomials P n (x) = M k (x; ao(r), c), where 

a -i\)\ (1.30) 
see jAkh65llKS98j . From this it follows that 



*V T^(fc E M fc (n; aa(r), c) |% for a" < 1 

k=0 



\n;a,a) r := < 



cr™ \n) r for a = 1 

(1.31) 

and the corresponding eigenvalues are 

:=2n + a (r). (1.32) 

Since the spectra of A and b ajCr (A ) coincide we define the unitary 
operator V a>a by 

U ,<t |«) r := |n;a,cr) r . (1.33) 

Using the fact that b ajCT (A_) and b 0jCT (A + ) satisfy the commutation 
relations (jl.2j) . we can check that they act on eigenvectors |n; a, a) r as 
weighted shift operators: 

b a , CT (A_) |n; a, cr) r = A/n(n + a (r) - 1) |n - 1; a, <r) r , (1.34) 
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b O)0 -(A + ) \n; a, a) r = yj(n + a (r))(n + 1) \n + 1; a, cr) r . (1.35) 

Thus the unitary operator U ajCr defined by (jl.33j) satisfies (jl.26|) on 
dense subset of 7i r consisiting of finite linear combinations of vectors 
\n\ b, ix) r for n G N U {0} and (6, 7r) G IR+ x Z 2 . In such way we have 
constructed the unitary representation of the group M+ x Z2, which 
implements the action ()1.25|) . However, it is not possible to extend it 
to whole 23 since b_ a cr can be decomposed into 

b-a,a = ba,<7 ° b-1,1, (1.36) 

where 

b_i,i(X) = -X (1.37) 

for XeA 

1.3. Integrable one-mode Hamiltonians. 

The aim of this section is to integrate quantum system described by 
arbitrary self-adjoint operator belonging to the multiboson algebra A: 

H Mi , := ^A + ^(A_ + A + ), (1.38) 

where (jm,u) Gl 2 \ {(0,0)}. 
Let us observe that 

\k) r = b k ^ \k - l) r + a k \k) r + b k \k+l) r , (1.39) 

where 

a k = ^^(2k + a (r)) (1.40a) 

b k = ^^y/(k + a (r))(k + l). (1.40b) 

If [i z/ the formula (jl.39j) is directly related to three term recurrence 
relation ()1.29|) for some family of orthonormal polynomials {P n } 
Since ^ ^ is divergent then there exists the unique measure duo on 
such the map F given by 

H r 3 \k) r i— > F(\k) r ) := P k G L 2 (R, duo) (1.41) 

is the isomorphism of Hilbert spaces with the property that 

F o H HHr o F" 1 = x, (1.42) 

where x is the operator of multiplication in L 2 (IR, duS), see |Chi78| 
OHT01 . Thus we gather that the spectrum spec(H A1!/ ) is the sup- 
port of measure duo. It means that by finding the measure duo and 



oo 

n=0- 



polynomials P n we obtain the evolution flow 

R 3 t — ► e itu ^ = F' 1 o e it& o F e Aut(H r ) (1.43) 

of quantum system described by the Hamiltonian H^ u . 

From definition of Bogoliubov group 03 it follows that transforma- 
tions b tt)Cr preserve the family of operators H^ u and the labels (/i, v) 
transform as follows 

(a, 1) : (fi,v) i — ^ (a [i,ais), ( a ? ~1) : (a*j ^) l— * ( az/ > a~V)- (1-44) 

This defines the action of the group 23 in the set M 2 \ {(0, 0)} of labels. 
We conclude that orbits of 03 are pairs of hiperbolae indexed by one 
real parameter eel 

O c := 23 • (c, 1) = {(X, y) G M 2 \ {(0, 0)} | xy = c}. (1.45) 

We can restrict our considerations to each component 7i r of decom- 
position ()1.13j) separately since they are invariant under the action of 
Due to the implementation formula (jl.26j) it is sufficient to find 
spectral decomposition for one operator from each orbit ([1.45)1 . e.g. 
H^ iV /h, tt^z-^z an d H 10 . Taking into account scaling by constant we 
can further restrict ourselves to three Hamiltonians H ljl; Hx,-i and 
Hi,o- 

Since Hi i = A then the spectral problem is trivial and is solved by 
formula (|1.15|) . 

For Hi o formula (jl.29|) defines Laguerre orthonormal polynomials 

P n {x) = L^~ l \2x) (1.46) 

and measure duj is given by 

dw(x) = 2{2x) ao{r) - 1 e' 2x 6{x)dx, (1.47) 

where 9 is Heaviside function. Thus the spectrum spec(Hi o) = M+ U 
{0}. 

For Hi ; _i formula (jl.29|) defines Meixner-Pollaczek orthonormal poly- 
nomials 

P n (x) = Pi~ } ( £ I ) (1.48) 



.2' 2 

and measure is du is given by 



du(x) = - 



1" ( -a (r) + i| 



2 

dx (1.49) 



and the spectrum spec(Hi i _i) = R, see |KS98j . 

In Appendix[X]we present complete list (for all /i, v) of Hilbert spaces 
L 2 (R,du) in which H MJy act as x. In Figure we have illustrated the 
decomposition of (/i, z/)-plane into sectors corresponding to different 
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families of orthogonal polynomials. Let us remark that in the first 
quadrant the spectra of H M!/ are discrete and bounded from below. On 
the boundary of this quadrant the spectra are R + U {0}. In the third 
quadrant all spectra are also discrete but bounded from above. On the 
boundary of this quadrant the spectra are M_ U {0}. In the second and 
fourth quadrants the spectra are R. 




Figure 1. Orthogonal polynomials assigned to H t 



1.4. Coherent state representation. 

The goal of this section is to give the coherent state representation 
of the Hamiltonians and flows e %m ^ generated by them. 

We will consider coherent states as eigenstates of A_, see 
Due to the decomposition ()1.13|) into irreducible representations, it is 
sufficient to restrict our considerations to each 7i r separately 

A-|0 r = ClC> r - (1-50) 

From ()1.16j) we see that the coherent states \C) r G H r are given by the 
series 

K),:=E 7= m \k) r (1-51) 
which converges for any ( 6 C and belongs to domain T>\. 
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The notion of the coherent states allows us to construct the anti- 
unitary embedding 

H r 3 |V>) i— > IrW(C) := (MOr & L 2 0(C,dfi r ) (1.52) 

of 7i r into the Hilbert space L 2 0(C, dp, r ) of holomorphic functions on 
C, which are square integrable with respect to the measure 

dUria) ■= Cr^if^ d P d ^ (1-53) 
27rr(a (^)) 

where ( = pe 1 ^ and K ao ( r ) is the modified Bessel function of the second 
kind, see |US97j . The space L 2 0(C, dp r ) has the reproducing kernel 



mO--=(v\0= ° F iL ( r ) (1-54) 
i.e. for any / G L 2 0(C, dfi r ) one has 

K(rj,Of(v)dfJr{v,v) = f{0- (1-55) 



The isomorphism I r gives the realization of the operators A , A + , A_ 
as the differential operators 

IroAooI^ 1 = 2C^ + a (r), (1.56) 
/,-oA+o/^ 1 = C, (1-57) 
I,oA_o/- = («d(r) + c|)| (1-58) 

acting in L 2 0(C, d/i r ). In order to describe them as the generators of 
the discrete series a (r) = 2, 3, . . . representation of the group SL{2, IR), 
let us consider a unitary integral transform 

V : L 2 0(C, dn r ) — ► L 2 £>(D>, du r ), (1.59) 

where © := {z E C | \z\ < 1} and 

dv r [z, z) := (1 — p| J a 2, (1.60) 

given by 

Vf(z) := | e^f(C)dfi r (z,z). (1.61) 

c 



11 

Using (jl.filj) we find that in the space L 2 C(D, dv r ) the operators 
(ll.56j) - (jl.58|) are given by 

V o I r .o A o r 1 oV' 1 = 2z^- + a (r) 

dz 

Po/.oA+o/^or 1 = z 2 ^- + a (r)z (1.62) 

dz 

v o i r o a o r 1 o v~ l = — 

az 

and they are the generators of the discrete series representation 

U?>W<j>(z) = (bz + a)- a «^<P (t^^ (1-63) 

of the group SL(2, R) in L 2 0(U),dv r ), see |Pe r86j . Here we have iden- 
tified SX(2,R) with SU(1, 1) using the isomorphism 

(1.64) 

From (jl.63j) we obtain the representation of SL(2, R) in L 2 0(C, dfi r ) 
in the integral form 

T; o(r) /(C) := | iC (r) (C,^) /fo) <W»7^) (1-65) 

c 

with the kernel 



In order to obtain the explicit formula for the evolution operators 
e iH M „t us consider the corresponding family of one-parameter sub- 
groups of SL(2,R) 



The infinitesimal generator of the action T^fl on L 2 0(C, d/i r ) is 
J r o (zH^t) o J,. -1 . Thus the time evolution operator I r o ( e * H M^*) J^. -1 
acts in L 2 0(C, dp, r ) as the integral operator (jl.65|) with the kernel 

C,W^'^ := / e v(i) ^^ (t) (-M*)™ + <V(*)J oV r (w,I(J), 

ID 

(1.68) 
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where 



/ 1 \ / r~ ~ i\ I isin(^7t) + 
<V,(£) = cos (v/i^ t) H 9 / (1.69) 



r /,n ^sin (^t) (-/i + i/) 

W*) = ^= • (L70) 



Ending let us remark that in L 2 0(C, dn r ) the Hamiltonian H M „ is 
represented second order differential operator 

- ^ H + + (■*> + 4) I) 

(1.71) 

and in L 2 (9(©, dv r ) as a first order differential operator 

PoJ r oH / ^o/- 1 op- 1 = ^ (^ z ± + a (r))+^ (V + 1)1 + a„(r)^ . 

(1.72) 

In Section 11.31 it was shown that this operator has discrete spectrum 
for \xv > The eigenproblem 

/ r oH F o/;V = A^ (1.73) 

in the Hilbert space L 2 0(C, d[i r ) leads to the degenerate Gauss equa- 
tion and its solutions can be expressed in terms of confluent hypergeo- 
metric function and Laguerre-type function, see |Smi61j . On the other 
hand the eigenproblem 

VoI r on flu oI; 1 oV- 1 cf) = \(f) (1.74) 

in the Hilbert space L 2 0(1D), dv r ) leads to the first order differential 
equation and its solution in the case /i > u > for the eigenvalue A is 
given by 

, (z) ,c(,- + ^^) A ( z+ Ut±^l) B , (1 . 75) 

where 



/i — v J \ /j, — v 



= - v) _ a (r) B = -\(n-u) _ a {r) 

2JJw 2 IJJw 2 1 ' ' 



Thus 

i) = (VcI r )- l cj> (1.77) 

and this kind of relationship gives the various integral representations 
of corresponding special functions. However the investigation of these 
problems is not the goal of the paper. 
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2. Two-mode systems 

2.1. Direct sum of multiboson sl(2, R) algebras. 

In this section we will work within two-mode bosonic Fock space 
7i£g>7i with the standard annihilation, creation operators a^, a*, i — 0, 1, 
acting in i th mode and with the orthonormal basis {|rio, ni =o 
consisting of eigenvectors of rij = a*aj. 

Let us consider the two-mode analogue of the investigation presented 
in the previous section. Let A © B = so{2,2) be the direct sum of 
two multiboson s/(2, E) operator algebras A := span{A , A_, A + } and 
B := span{B , B_, B + } acting on V <g> V Q C 7i ® H, where 

A := a (n ,ni) A := a_(n ,ni) a '° A + := (a*) /o a_(n , m) 

(2.1) 



B := Aj(no,ni) B := /?_(n , ni ) a/ 1 B + := (a*)^_(n , m) 

(2.2) 

and A and B satisfy (OJ) and (Ol) 

Proceeding as in previous section we gather that the commutation 
relations for A®B imply the difference equations on functions ao, 
/3o, (3-. The solution of these equation gives the following expressions 
on operators: 

A = T (n - Ro) + a (Ro, Ri) (2.3a) 



n 



Q-(n - R ) + a (Ro, Ri)) f^(n - R ) + a 



lo 



B 



;2.3b) 



B = y (nx - Ri) + /9 (Ro, Ri) (2.4a) 



g (ni - R0 + AtRo.R,)) (i(n, - Rl ) + l) 



(2.4b) 

where Ri, i = 0,1, are defined by ()1.9|) for i th mode and cto, A) are arbi- 
trary positive functions on the product of spectra (spec R ) x (spec Ri). 

We observe that Rj commute with A®B and that Casimir operators 
Qaj (see ()1.20j) ) can be expressed as functions (not invertible in 
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general) of Ro and Ri. Thus we will decompose two-mode Fock space 

lo-i h-i 

ro=0 ri=0 

into common eigenspaces of R« 

H ro , ri ■= span{ \k , fa) ro n := \k l + r , Mi + n) \ k , fa G N U {0}}, 

(2.6) 

on which A® B\n TOtr ac ^s irreducibly. 

Comparing the formulae (|2.Hjl - (|2.4jl to (jl. 12)1 we extend all generators 
of A © B to the common domain T>i ® T>\. 

2.2. Hamiltonians. 

Let us consider a class of sl(2, M) subalgebras X> := span{D , D_, D + } C 
A © B defined by 



D 

D_ 



o 



b a , CT (Ao) + b feiT (B ), (2.7) 
b aj(T (A_) + b 6)T (B_), 
b a , (T (A+) + b 6 , T (B+), 



for (a, a), (6, r) e 03. For simplicity we have omitted indices in D , 
D_, D + corresponding to the elements of 03. 

Since D , D_, D + satisfy (|1.2|) . the Casimir for this subalgebra is of 
the form 

C v = iD 2 - D D + - D + D (2.8) 

and it commutes with the Casimirs C_4 and Qg. 

Main goal of this section is to integrate the two-mode bosonic quan- 
tum system described by the Hamiltonian 

H := ~ (C v -C A - C B ) . (2.9) 

Casimir operator Ct> is not defined on T>\ © T>\ but due to the explicit 
formula 

H - <^lA B - ^ ^ , 2 ,0) 

a 2 -b 2 (a + b) 2 

+ ^r(A.B- + A B + )-.rL^(A + B_ + A_B + ) 

we see that the symmetric operator H is well defined on T)\ ®D\. 

Besides obvious integral of motions Ro, Ri this Hamiltonian also 
commutes with Dq which is so called Manley-Rowe integral of motion 
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for H. It means that the investigation of the the spectral decomposition 
of H can be reduced to the eigenspaces of D in each TC rQ>ri . 

Due to the implementation formula ()1.26|) we can unitarily transform 
H by <g> id + id ® U| 6 |-t t to the Hamiltonian 

Hd := ^A B + A + B + + A B (2.11) 

for ab > and to the Hamiltonian 

He := ^AqBo + AB+ + A+B (2.12) 

for ab < which correspond to indices (a, cr),(b, r) in ()2.7j) being 
(1, -1), (1, 1) or (1, -1), (-1, 1) respectively. 

Spectral decompositions of self-adjoint extensions of (J2.11j) and (|2.12j) 
are essentially different thus we will discuss it separately. In order to 
simplify the notation, we will replace symbols ao(r ,ri), /3o( r (b r i) by 
«o and (3o in all formulae below when there will be no risk of confusion. 

2.3. Spectral decomposition of Hp. 

In this case the operator Do assumes the form 

D = A + B . (2.13) 

Thus its spectrum consists of numbers 2K + a + O for K = 0, 1, . . . 
and corresponding eigenspaces 

< /ri := span{ |fc)* n | k = 0, 1, . . . , K} (2.14) 

are (K + l)-dimentional, where 

\k)«, ri :=\k,K-k) ron . (2.15) 

So we get the decomposition 

oo 

^ro,n = ^^^ro.rr (2.16) 
K=0 

onto subspaces invariant with respect to Ho. Hamiltonian Hd| w k r ^ is 
self-adjoint and thus admits unique self-adjoint extension in H.®TL. 
Moreover it acts in basis {\k) ro ri } k=0 as follows 

lCn = h-i \k - 1>* n + a k \k)« !n + b k \k + l)£ n , (2.17) 

where 

a k = ^(2k + a )(2(K-k)+{3 ), (2.18) 

b k = y/(k + l)(k + a )(K -k)(K-k + fa). (2.19) 
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Thus, analogously to the Section II the problem of finding the spec- 
tral decomposition of can be solved by means of some family of 
orthogonal polynomials. In this case the eigenvalues of HdI^k are 
given by 

E° :=n(n + a + po-l) + ^a Po (2.20) 
for n = 0, . . . K and corresponding eigenvectors assume the form 

K 

iD\K 



ro,rx 

k=0 



where 



P,(„) := /?, ( - ±a (3 ;a - 1, A, - ) = (2.22) 



a - 1 + k\ fp + K - k\ f-k, -n, n + a + O - 1 
k \ K-k 32 \ a ,-K 



are normalized Dual Hahn polynomials, see |KS98j . The polynomials 
{Pk}k=o form the orthonormal basis in L 2 (M, duS) with the measure 
given by 



K 



xS I x - E® - ^a /?o ) cfe. 



n=0 

1 

2 

Finally we get the following spectral decomposition 

l -l h-l OO K , , „£>, 

h d = EEEE^§#/- < 2 - M > 

r =0ri=0if=0 fc=0 \ n\ nl 

In particular we see that the eigenspaces of mTi^Ti are infinite 
dimensional. 

2.4. Spectral decomposition of He- 

In this case the operator Do assumes the form 

D = A - B . (2.25) 

Thus its spectrum consists of numbers 2K + a — A) for K G Z and 
corresponding eigenspaces 

< iri :=span{|fc)^ ri |fc = 0,l,...} (2.26) 
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are infinite dimentional, where 

\k)r ,n := 
So we get the decomposition 



K + k,k) rori forfT>0 
k,k — K)'. for K < 



(2.27) 



(2.28) 



onto subspaces invariant with respect to H^. The Hamiltonian Hc| w k 



also acts in basis {\k)^ o ri }^L like in previous section, i.e.: 

He = h-i \k - l)£ n + a k |fc>* n + b k \k + l)£ n , (2.29) 

where now 



-i(2(K + k) + a )(2k + /3 ) for K > 
-|(2fc + ao)(2(A;-liO + A)) forir<0 



(2.30) 



V(AT + fc + a )(K +k+ l)(k + P )(k + 1) for^>0 
■y/(k + a )(k +l)(k-K + p )(k - K + 1) for AT < ' 

(2.31) 

Three-term recurrence relation (jl.27|) with these coefficients is solved 
by normalized Continuous Dual Hahn polynomials 

P n {x) = S n (-x; u, v, w) := (2.32) 

(u + V) n (u + W) n ^ ( —fl U + y/x U - 



xjT{n + u + v)T(n + u + w)T(n + v + w)n\ 3 2 V " - " + »' 
which are orthogonal with respect to the measure 



du(x) 



Y(u + i\J—x)T(v + i\J—x)T(w + i\ 



-x) 



r(2ii 



-X) 



1 



-x 



Y{u + v)Y{u + w)Y{v — u)Y{w — u) 



[—x)dx+ 
(2.33) 



n=0,l,2,.. 
u+n<0 



r(-2o) 

(2w) n (w + l) ra (tt + t>) n (n + w) n | 
(li)n(lt - v + l) n (u -w + l) n n\ 



X 



■l) n 8(x- (u + n) 2 )dx, 



see |KS98j . The parameters u, v, w depend on the constants a , j3 
and K in the following way. For K > we have 

±(/? -a + l) if A) - «o e (-oo, -1) 

|(ao + A)-l) if/? -a G (-1,2^ + 1) 

if+ §(a - A) + 1) if A) - «o e (2K+ l,oo) 
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K +1(00-/30 + 1) 
Upo-oo + l) 
Uao + po-1) 



\(ao + (3o-l) 
tr = < K + |(a - 0o + l 
i(/3 -Q!o + l) 
and for 7^ < we have 

-K +\{(3 -ao + l) if/9 -a e 



if ft-«o6 
if p - a E 
if (3q - a e 

if p - «o e 
if /3 - a e 
if p - a e 



u 



Uuo + po-l) 
\{ao-po + l) 



if 0o 

if A) 



|(a -/? + l) 
-K+K^o-ao + l) 
|(a + A)-l) 

|(a + A)-l) 
|(a -A) + l) 
-lf+i(A)-ao + l) 



a Q e 



if - "o G 

if /? - a G 
if /? - a E 

if /3 - a E 
if /3o - a G 
if /? - a E 



-oo, -1) 
-1,21T+1) 
21sT+ l,oo) 

-00, -1) 
-1,2K + 1) 
2K + l,oo) 

-oo,-l +2ZT) 
-1 + 2K,1) 
l,oo) 

-00,-1 + 210 
-1 + 2K,1) 
l,oo) 

-00, -1 + 210 
-1 + 2K,1) 
l,oo). 



Thus Hc| w k admits unique self-adjoint extension (see |AG93j ). so 
He also admits unique self-adjoint extension 'va.7i®7i. The operator 
F o H c \ h k r o F^ 1 , where F is defined like in (11.41)1 . acts in L 2 (IR, duo) 

as operator of multiplication by x — |((a — l) 2 + (Po — I) 2 — l) and 
therefore the spectrum of Hc\n K is 



spec HcI-wk 



1 



U 



n) 



—00, 



1} 



((oo - l) 2 + (A) - I) 2 - 1) U (2.34) 



(A) 



1) 



0,1, 



A u + n < 



so it always consists of the continuous part (unbounded from below) 
and if u < there are — [u] points in the discrete part. The eigenvectors 
corresponding to discrete part are of the form 



K 

ro,n 



00 

E 

fc=0 



Pfc^ + r*) 2 ) |fc) 



K 

ro,n 



(2.35) 



for n = 0, 1, . . . , — [u]. 

We see that also in this case the eigenspaces of H c in 7i ® 7i (if they 
exist) are infinite dimensional. 

Similarly to the one-mode case one has the coherent state represen- 
tation for the Lie algebra A © B and thus also for the Hamiltonian 
(12.9)) . This representation is obtained as the tensor product of the one- 
mode coherent state representations. Hence all formulae concerning 
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the two-mode coherent state representation are easily obtained by the 
tensoring procedure. Therefore we will not present it here. 

3. Physical remarks 

In the previous chapters we have found spectral decomposition of 
the Hamiltonians H Mi , (jl.38|) and H ()2.9|) . i.e. we have integrated the 
considered systems. 

The aim of this chapter is to discuss the physical interpretation of 
the Hamiltonians and H. First of all let us state that we will 
interpret these operators as Hamiltonians of one-mode and two-mode 
bosonic quantum fields in non-linear medium. 

We consider H (resp. H^) as interaction Hamiltonians for the sys- 
tem and, as shown in the paper HOT02J, we introduce full Hamiltonian 
by the formula 

H F := H + e" iHo 'H e iH °*, (3.1) 

where H := Wono+Wini (resp. H := un) is free bosonic Hamiltonian. 
The solutions of Schrodinger equation are given by 

[V(*)> = e- iHot e~ iHt \ifj(0)} (3.2) 

and allow us to compute the time evolution of the physical character- 
istics of the system, e.g. mean numbers of particles, their dispersion, 
Fano factors, correlation functions, and squeezing factors. 

To give more explicit interpretation of H let us observe that from 
the formulae (l2~TU|) and (JO), it follows that H is of the form: 

H = £oo(n , ni) + g+-(n , ni )(a*) ,0 a£ + g_ (n , n 1 )4°+ (3.3) 

+flr -(n ,ni)a{ 1 + ( n o, n^a^af 1 + /i.e.. 
The term goo( n cb n i) corresponds to Kerr- type effects in the non-linear 
medium (e.g. bistability), see |PL98j . The term <? + _(n , ni)(a£) '°a 1 1 
describes intensity dependent (parametric) conversion of the cluster of 
lo bosons in the first mode into the cluster of l\ bosons in the second 
mode, i.e. the factor (aj) z °a^ 1 describes the absorpion by the medium 
of one bosonic cluster and emission of another cluster, while the factor 
(7 + _(n ,ni) is intensity dependent coupling constant. The conjugated 
term describes the reverse conversion of clusters. By analogy the re- 
maining terms describe the process of intensity dependent absorption 
of the cluster of bosons in one or both modes. The conjugated terms 
describe the emission of the same clusters. The process of emission- 
absorption occurs with the probability depending on the functions gn, 
where k, I = 0, +, — , so it essentially depends on the number of bosons 
in the medium. 
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For ab > the spectrum of H is purely discrete thus the considered 
system has only bound states. On the other hand for ab < the 
spectrum of H has continuous part and (in particular cases) discrete 
part. Thus the considered system has scattering states and (in some 
cases) bound states. Let us remark that due to the fact that in this 
case the spectrum is unbounded from below, it is correct from physical 
point of view to consider Hamiltonian — H. 

Some representatives of integrated class of Hamiltonians, related to 
Examples 11.11 and 11.21 are: 

H, = n + n x + 2n + a 2 a 2 + (a*) 2 (a*) 2 , (3.4) 
H// = n + m + 2n ri! + a 2 (a*) 2 + a 2 ^) 2 , (3.5) 
H//7 = n + ni + 2n ni + y/n^ a* a? + Vn + 1 a (a^) 2 , (3.6) 

Hjy = no+ni+2110 111 + v/noni a^aj + v/ (n + l)(ni + 1) a a!. (3.7) 

Detailed physical analysis of systems described by Hamiltonians Hm 
and H/v is the subject of the paper jTOH+06j . 

Physical interpretation of Hamiltonians H Mi , is analogous, i.e. they 
describe parametric absorption-emission of one-mode of bosons in non- 
linear medium. 

Concluding we see that investigated Hamiltonians include a wide 
range of non-linear processes of interaction of bosons with medium. In 
particular applications they can describe the models of physical systems 
consisting e.g. of photons, phonons, magnetons, or Cooper pairs. 

Appendix A. Spectral decomposition of one-mode 

Hamiltonians 

In this Appendix we present formulae for spectral decomposition of 
Hamiltonians (jl.38|) . The problem splits into several cases. 

Case 1. v = 0, fi ^ 

H Mi , is related to Laguerre orthonormal polynomials P n (x) = Ln ( 
In this case 

and spectrum spec(H /ii/ ) = M + U {0} for /x > and spec(H AtJy ) = 
R- U {0} for n < 0. 
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Case 2. /i = 0, v ^ 

H Ml/ is also related to Laguerre orthonormal polynomials but by dif- 
ferent formula P n (x) = (-l) n Li ao(rM) (f ). In this case 



duj(x) — — ( — 



2 /2x\ ao(r) - 1 2 , 



e-~6 -)dx 



(A.2) 



and spectrum spec(H^) = R + U {0} for v > and spec(H^) 
R_ U {0} for // < 0. 



Case 3. \x > 0, v < 

H MZ , is related to Meixner-Pollaczek orthonormal polynomials P n {x) 



Pn 



==; 0), for = arccos(— ^J). In this case 



dcu(x) 



(20-tt) 



and spectrum spec(H Mi ,) 



(A.3) 



Case 4. ^ < 0, z/ > 

H MZ , is related to Meixner-Pollaczek orthonormal polynomials P n (x) 

P«. 2 L r= ; 0), for = arccos(— jfr^)- In this case 
I 



duj(x) 



(20-7r)-=2= 



and spectrum spec(H Mi ,) = M. 



r (r o(r) - 2 v^ 



<ix (A.4) 



Case 5. fi > v > 

is related to Meixner orthonormal polynomials P n (x) = M n (j-^=- 



ao(r) . 
2 



;a (r),c), for c 



In this case 



(icj(x) = 5(x — a G {r)^/Jw — n) - °^^ n c n d x 



n=0 



(A.5) 



and spectrum spec(H Mi ,) = n + a ( r )\//^' I n = 0,1,2,.. .}. 

Moreover the set of vectors 



(a (r)) n c 



-J2M k (n;a (r),c) \k\ 



(A.6) 



fc=0 



n=0 



is an orthonormal basis in H r consisting of eigenvectors of H M „. 
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Case 6. /i < v < 

is related to Meixner orthonormal polynomials P n (x) = M n (^=- 

2^1; a„(r), c), for c = In this case 

du(x) = ^5{x + a (r) v T^ + n) ^°^^ n c n cfe (A.7) 

n=0 n ' 

and spectrum spec(H M1/ ) = {— l^fjw n — ao{r)^/Jw \ n — 0, 1, 2, . . .}. 
Moreover the set of vectors 



nl 



(ao(r))nC nZ M ^Mr),c) \k)A (A.8) 

A;— ) n=0 

is an orthonormal basis in Ti r consisting of eigenvectors of H^. 

Case 7. v > /i > 

is related to Meixner orthonormal polynomials P n (x) = (— l) n M n ( 

a (r), c), for c = g^/jE - In this case 

dw(a;) = ^ <J(x - a (r)v^ - 2^ n) ^°^ w c w cto (A.9) 

n=0 

and spectrum spec(H Mi ,) = n + ao( r )v^ 7 I n = 0,1,2,.. .}. 

Moreover the set of vectors 



nl 



J2(-l) k M k (n;a (r),c) \k) r \ (A.10) 



(an(r)) n c n 

is an orthonormal basis in H r consisting of eigenvectors of H^. 

Case 8. v < \x < 

is related to Meixner orthonormal polynomials P n (x) = (—l) n M n 

a (r), c), for c = ^ +2 /E . I n this case 

du(x) = ^J(x + a (r) v ^ + 2 v ^ir;n)^^c n rfx (A.ll) 

n=0 

and spectrum spec(H^) = {— l^fjw n — ao{r)^fJw \ n — 0, 1, 2, . . .}. 
Moreover the set of vectors 



f - Mj _g(-l)*M t ( n;aoW ,c) |*>,j^ (A.12) 
is an orthonormal basis in 7i r consisting of eigenvectors of H^. 
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Case 9. \l = v 

In this case formula (jl.39j) means that H MJ , is diagonal in Fock basis 
and spectrum spec(H^) = {2/in + fj.a (r) \ n = 0, 1, 2 . . .}. 
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